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In examining a select sample of over 200 blazars of known jet kinetic power {L]^in) ctnd well-characterized SEDs, 
we found 1] that Intermediate synchrotron-peaking (ISP) blazars may have lower gamma-ray output than high 
synchrotron-peaking (HSP) blazars of similar Lkin, consistent with our hypothesis that ISP blazars are less- 
beamed versions of HSP blazars, rather than a distinct population. Further, by using the radio core dominance 
as a measure of relative beaming, we find that gamma-ray luminosity depends on beaming in a consistent way 
for blazars ranging over all jet kinetic powers (10^^ — 10^^ ergs s"-*^). We re-examine the gamma-ray properties 
of this core sample of blazars using the 1-year LAT catalog We find that for weak jets, the ratio of inverse 
Compton to synchrotron emission remains constant with increased beaming, consistent with an SSC model for 
the jet emission, while the most powerful jets show a strong increase in Compton dominance with orientation, 
consistent with an external Compton (EC) emission model. 
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1. Introduction 

Radio-loud active galactic nuclei (RL AGN) are be- 
lieved to contain an actively accreting super-massive 
black hole (10^ - 10^^ ^0) which generates bi-polar 
jets of relativistic material reaching up to Mpc in 
scale and with luminosities up to 10^^ ergs s~^, pow- 
erful enough to heat the intra-cluster medium and 
potentially produce a feedback mechanism in galaxy 
formation [3]. When these relativistic jets are seen 
with the axis near to the line of sight, they are 
called blazars (see [4] for a review). The jet radia- 
tion is strongly Doppler-boosted and dominates over 
the emission from the host galaxy, the dusty torus 
(thought to surround the nuclear black hole), the ac- 
cretion disk, and, except at very low frequencies, the 
isotropic emission from the giant lobes of gas fed by 
the jet over time. Thus blazars are naturally interest- 
ing objects as the magnified view of the jet emission 
can allow us to probe the structure and nature of these 
jets and ultimately to understand the physical mech- 
anisms behind their origin. 

The broad-band, highly luminous and variable 
spectrum from blazars is usually characterized as a 
'double-peaked' one, with one broad emission com- 
ponent believed to be synchrotron radiation peaking 
from sub-infrared energies to X-rays, and a second 
peak from inverse Compton emission at > MeV ener- 
gies. This second high-energy peak in blazars can be 
explained with both synchrotron-self Compton (SSC, 
in which the particles in the relativistic jet up-scatter 
synchrotron photons, [5|, |6|) and external Compton 
(EC) models, in the latter case with photons from 
either the accretion disk, broad line region (BLR), 
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or molecular torus 0, H, 0). Generally low -power 
and lineless objects are believed to radiate by SSC at 
high energies, while individually, many high-powered 
FSRQ have been more satisfactorily fit with EC mod- 
els [lo|, El. However a consistent basis for which 
sources require EC has not been demonstrated for any 
particular class of blazars. 

The location of the GeV emission is also a matter of 
active debate, as it comes from locations close to the 
central engine that remain unresolved, even with very 
long baseline interferometry (VLBI). If the emission 
region is located within the sub-pc scale BLR, the 
GeV emission of blazars with strong lines is due to 
EC scattering of BLR photons [12], which can lead to 
more efficient cooling and thus lower peak frequencies, 
as seen in the blazar sequence [13]. However, others 
agree that the high-energy emission may be coming 
several pc downstream as suggested by some multi- 
wavelength observations of variability [13, IHI. The 
EC versus SSC origin of the high-energy emission is 
an important question as it can be used as a diagnos- 
tic for the location of the gamma-ray emission, and 
the inferred structure of the jet. 

The Fermi gamma-ray telescope is giving us an un- 
precedented look at the high-energy emission from 
blazars, with nearly 700 blazar associations in the 1- 
year catalog (1-LAT)[2]. With such a large sample, we 
can begin to put some order in the phenomenology of 
blazars at high energies. In Section [2] we introduce the 
'Blazar Envelope' as seen by Fermi and discuss the im- 
portance of the jet kinetic power and orientation of the 
jet in determining the gamma-ray luminosity. In Sec- 
tion [3] we discuss a recent finding suggesting that the 
high-energy emission of high-power jets is dominated 
by EC rather than SSC. In Section 21 we summarize 
our findings. 
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Figure 1: The synchrotron Lpeak - J^peak plane, in which 
the blazar sequence first found an ant i- correlation 
between bolometric luminosity (with Lpeak) and the 
synchrotron peak frequency. Using only well-sampled 
SEDs reveals a possible alternative to the continuous 
sequence, in which two populations show dramatically 
different behaviors in this plane. Track (A) is the 
de-beaming trail in synchrotron peak for a standard 
single Lorentz- factor jet ('strong' type), track (B) shows 
the markedly horizontal movement typical of a 
decelerating jet model characterized by velocity gradients 
('weak' type). Adapted from J.]. 
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Figure 2: The Inverse Compton Envelope. Figure is 
adapted from [16]. The Fermi-detected AGN (large, dark 
blue circles and magenta squares) show that alignment 
plays a strong effect in the gamma-ray output. 
Interestingly, while the FR II sources appear to drop 
directly below the powerful FSRQ sources, FR I radio 
galaxies appear to follow a more horizontal track, similar 
to our findings in the synchrotron envelope (Figure 1, 
left). The presence of a forbidden zone (empty region at 
upper right) suggests that there is a sequence in the IC 
plane, similar to that found in the synchrotron (Figure 



2. The Inverse Compton Envelope as 
seen by Fermi 

The relativist ic jets in RL AGN are generally con- 
fined to a small opening angle and thus most blazars 
are seen within a small range of orientation angles 
from the jet axis. It is often assumed that bright 
blazars in flux-limited samples are "well- aligned" and 
therefore the exact viewing angle is of little impor- 
tance (e.g., for modeling, often the assumption = 
l/F is used). However, as surveys go deeper, it is 
expected that more and more relatively misaligned 
sources will be observed. 

Fossati et al. (1998) demonstrated an anti- 
correlation between the synchrotron peak luminosity 
and the peak frequency, forming the now canonical 
'blazar sequence' of (presumably) well- aligned sources 
[13 • Lower-luminosity sources appearing in the space 
below the blazar sequence are expected due to some 
sources being less aligned to the line of sight. Indeed, 
[l8[ found that new sources they identified modify the 
blazar sequence to an envelope, with the area below 
the blazar sequence populated with sources. Similar 
envelopes were found by Anton & Browne and Niep- 
pola et al. [iH, [20|. In recent work, we filled the 
synchrotron Lpeak — ^peak plane using a much larger 
sample of jets over a wide range of orientations (from 
blazars to sources misaligned enough to be seen as ra- 
dio galaxies), with much-improved SED sampling f]\. 



This was the first work to consider both the jet ki- 
netic power (Lkin, as measured from the isotropic ra- 
dio emission) and radio core dominance (i?, a measure 
of jet orientation), as important factors. Strikingly, we 
found that the blazar sequence is actually broken into 
two populations in the Lpeak — ^peak plane (Figure 
[1]). The 'weak' jets consist entirely of sources with 
Lkin <10^^"^ ergs s~^, and extend from mis- aligned 
FR I radio galaxies at low Vpeak up to well-aligned 
(as measured by R) HSP BL Lacs. The 'strong' jets 
comprised a population of low i^peak jets which drop 
in peak luminosity rapidly with decreasing R. The 
critical transition between weak and strong jets ap- 
pears to occur at an (dimensionless) accretion rate of 
rhcr ^ 10~^ (see also [21, 22]), as estimated from the 
ratio Lkin/LEdd^7 matching the divide suggested by 

We note that the large number of blazars seen by 
Fermi, as well as the detection of several radio galax- 
ies [16], implies that we should have a large range in 
orientations in the Fermi samples. As see in Figure O 
an 'envelope' is indeed seen in the plane of gamma-ray 



-•^The Eddington luminosity was calculated using the average 
black hole mass estimate from the literature, as will be discussed 
in a forthcoming publication. 
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Figure 3: Total gamma ray output versus radio core 
dominance for the 130+ sources with good SED coverage 
detected by Fermi in the 1-LAT. Jet power has a very 
strong effect on the gamma-ray luminosity, but when 
sources are sorted by their jet powers, the dependence on 
radio core dominance also becomes apparent. The slopes 
shown are statistically significant at the 95% level. 



luminosity (> 100 MeV) versus the LAT photon in- 
dex (which approximately tracks the peak frequency) . 
There is clearly a 'forbidden zone' at the upper right, 
with the edge of the envelope presumably comprised 
of the most aligned objects. The sources are divided 
into the typical higher-power and broad-lined flat- 
spectrum radio quasars (FSRQ) and the lower-power, 
lineless BL Lac objects (BLL). While this figure lacks 
the precise information needed to find the two popu- 
lations suggested by [1], it is clear the the FR I radio 
galaxies (which are unified with most BL Lacs) appear 
at much lower IC peak frequencies, forcing a similar 
horizontal track from the low-power, aligned BLLs, 
just as was found in the synchrotron plane. It is in- 
teresting to note that the FR II radio galaxies, which 
would be generally unified with the FSRQ population, 
are both directly beneath them in peak frequency, but 
are also not particularly mis-aligned, which suggests 
that these sources might drop out of Fermi detection 
more quickly than the weak sources. 

It is interesting to consider the effect of jet kinetic 
power as well as orientation on the total LAT-band 
luminosity. As shown in Figure [3l as a total popula- 
tion the IFGL blazars appear as a scatter-plot in a 
plot of versus R. However, when the jet kinetic 
power is shown (color scale), it becomes clear the 
increases with both radio core dominance (i.e., angle 
to the line-of-sight) and, even more dramatically, with 
the jet kinetic power. We note that Lkin is calculated 
from the low-frequency isotropic emission and is there- 
fore a completely independent measurement from L^, 
making it unlikely for this trend to be the result of 
the selection effects of the 1-LAT sample. 



3. Evidence for External Compton in 
High-Power Jets 

With a large sample of blazars spanning a range 
of orientations as shown in the previous section, it is 
possible to look for signatures of different high-energy 
emission mechanisms in a population of blazars. The 
Doppler beaming factor, (5, is a function of the Lorentz 
factor r and the orientation angle, and the effect of S 
on the observed (monochromatic) luminosity is: 

L = Lo<53+" (1) 

where Lq is the rest-frame luminosity (at 5 = 1), 
and a is the energy spectral index at the frequency of 
interest. The exponent 3-\-a is the value assumed for 
a 'moving blob' in the jet. If the emission comes from 
a standing shock, the exponent would be 2+a [24]. 

In the case of emission by SSC, the IC peak has a 
beaming pattern which is identical to the synchrotron 
(i.e., it follows Equation [1]) . However, for EC models, 
the beaming at high energies goes as Lie = Lo,ic^^^^^" 
[25|, |26|. The larger exponent indicates that as a 
source with significant EC emission is aligned, the IC 
peak should be more and more dominant over the syn- 
chrotron peak (i.e., the ratio Ljc/i^sync will increase 
with R. We can measure the Compton Dominance 

RCD = log(Lic/Lsync) Oclog (^(3+«)/(4+2c.) ) ^ 

that the dependence on S goes with an exponent of 2 
when comparing the IC and synchrotron peaks (where 
a = 1). At radio frequencies, where the spectral in- 
dex is confined to a fairly narrow range in values of 
ar = — 0.5, the beaming exponent will be ~ 3 - 3.5. 

In Figure [H we examine the relationship between 
RcB and the radio core dominance, i?, for two sub- 
samples of the 1-LAT blazars. The pink triangles are 
high-jet-power sources (log Lkin > 45.5), while the 
blue are low-power (log Lkin < 43.5). The flat dis- 
tribution in the latter group is consistent with the 
SSC models generally used to model these types of 
sources. For the high-power group, however, we find 
a significant trend of increasing Compton dominance 
with increasing radio core dominance (i.e., alignment). 
This can be interpreted as a signature of an EC process 
for a population of jets characterized by high kinetic 
powers. 

However, from the above relations, the expected 
slope of the correlation should be close to 2/3, up 
to a max of 1 under reasonable assumptions. The 
slope in Figure |4] is clearly much higher, in fact con- 
sistent with a value of ~ 2. While a more detailed in- 
vestigation is in progress, we note that this high slope 
might mean that the assumption that the Lorentz fac- 
tor of the plasma emitting in the gamma-rays is the 
same as that emitting at GHz frequencies is incor- 
rect. If a roughly 2:1 ratio is assumed (i.e., Tjc = 15 
while Fradio = 8), the higher slope can be explained. 
This is consistent to some extent with the difference in 
variability timescales for very high energies (minutes 
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Figure 4: We compare strong- jet and weak- jet sources in 
terms of their inverse Compton dominance over 
synchrotron emission. We find that for the strong jets, 
the inverse Compton dominance increases dramaticahy 
(up to a factor of 100) with increasing ahgnment, as 
measured by radio core dominance. Conversely, the weak 
jets show no such increase. The IC peak is measured by 
fitting the X-ray, gamma-ray, and TeV (where available) 
points with a parabola and taking the peak luminosity. 
(Dashed line shows slope=2 for reference). 



to hours) versus radio (much longer), which indicates 
that the radio emission is coming from a region which 
is slower and/or at larger scales. 

4. Conclusions 

We have shown that the Fermi gamma-ray satellite 
sees relativist ic jets in RL AGN over a wide range 
of orientation angles, including sources seen as ra- 
dio galaxies. There is a 'forbidden zone' of high- 
power, high-peak sources which remains empty, con- 
sistent with a spectral sequence similar to that seen 
in the synchrotron peak luminosity - peak frequency 
plane. Further, the much lower IC peak frequencies 
found for the few detected FR I radio galaxies sug- 
gest that the low-power class of FR I/BL Lacs de- 
beam along more horizontal tracks in the synchrotron 
and IC planes (Figures 1 and 2). We show that the 
total gamma-ray band luminosity depends on both 
orientation (measured through radio core dominance) 
and more heavily on the jet kinetic power (as mea- 
sured from low- frequency, isotropic radio emission). 
Finally, we present the first collective evidence for EC 
process in a sample of high-power jets (log Lkin > 
45.5), showing that as alignment (radio core domi- 
nance) increases, the Compton dominance increases 
dramatically. 

Acknowledgments 

This research has made use of the SIM- 
BAD database, operated at CDS, Strasbourg, 



France ( |http: / / simbad.u-strasbg.fr / simbad/ ) , and the 
NASA/IPAC Extragalactic Database (NED) which 
is operated by the Jet Propulsion Laboratory, Cal- 
ifornia Institute of Technology, under contract with 
the National Aeronautics and Space Administra- 
tion (http://nedwww.ipac.caltech.edu/). We utilized 
package np in the R language and environment for sta- 
tistical computing. GF and EM acknowledge support 
from NASA grants NNG05GJ10G, NNX06AE92G, 
and NNX09AR04G, as well as SAO grants G03- 
4147X and G05-6115X. MG acknowledges support 
from the NASA ATFP grant NNX08AG77G and 
NASA FERMI grant NNH08ZDA001N. The MO- 
JAVE project is supported under National Science 
Foundation grant 0807860- AST. 

References 

[1] Meyer, E. T., Fossati, G., Georganopoulos, M., & Lis- 
ter, M. L. 2011, ApJ, 740, 98 

[2] Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010, 
VizieR Online Data Catalog, 2188, 80405 

[3] McNamara, B.R. and Nulsen, P.E.J. , 2007 Ann. Rev. 
Astron. & Astrophys. 45, 117 

[4] Urry, C. M., & Padovani, P. 1995, PASP, 107, 803 

[5] Maraschi, L., Ghisellini, G., k Celotti, A. 1992, ApJ 
Lett., 397, L5 

[6] Marscher, A. P., & Travis, J. P 1996, A&As, 120, 537 
[7] Ghisellini, G., & Madau, P. 1996, MNRAS, 280, 67 
[8] Blazejowski, M., Sikora, M., Moderski, R., & Made- 

jski, G. M. 2000, ApJ, 545, 107 
[9] Sikora, M., Stawarz, L., Moderski, R., Nalewajko, K., 
& Madejski, G. M. 2009, ApJ, 704, 38 
Ghisellini, G., Tavecchio, F., Foschini, L., et al. 2010, 
MNRAS, 402, 497 

Vercellone, S., Striani, E., Vittorini, V., et al. 2011, 
ApJ Lett., 736, L38 

Sikora, M., Begelman, M. C, & Rees, M. J. 1994, 
ApJ, 421, 153 

Ghisellini, G., Celotti, A., Fossati, G., Maraschi, L., 
& Comastri, A. 1998, MNRAS, 301, 451 
Bottcher, M., Reimer, A., & Marscher, A. P. 2009, 
ApJ, 703, 1168 

Agudo, L, Marscher, A. P., Jorstad, S. G., et al. 2011, 
ApJ Lett., 735, LIO 

Abdo, A. A., Ackermann, M., Ajello, M., et al. 2010, 
ApJ, 720, 912 

Fossati, G., Maraschi, L., Celotti, A., Comastri, A., 



[lo: 
[11 
[12: 
[13: 

[14 

[is: 
[16: 

[17: 
[is: 

[19 

[20: 

[21 

[22' 
[23^ 

[24 
[25^ 
[26^ 



& Ghisellini, G. 1998, MNRAS, 299, 433 

Padovani, P., Perlman, E. S., Landt, H., Giommi, P., 

& Perri, M. 2003, ApJ, 588, 128 

Anton, S., & Browne, 1. W. A. 2005, MNRAS, 356, 
225 

Nieppola, E., Tornikoski, M., & Valtaoja, E. 2006, 
A&A, 445, 441 

Georganopoulos, M., Meyer, E. T., Fossati, G., & Lis- 
ter, M. L. 2011, arXiv:1111.4711 
Ghisellini, G., & Celotti, A. 2001, A&A, 379, LI 
Narayan, R., Garcia, M. R., k McClintock, J. E. 1997, 
ApJ Lett., 478, L79 

Lind, K. R., k Blandford, R. D. 1985, ApJ, 295, 358 
Dermer, C. D. 1995, ApJ Lett., 446, L63 
Georganopoulos, M., Kirk, J. G., k Mastichiadis, A. 
2001, ApJ, 561, 111 



eConf C110509 



